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ABSTRACT 

We present the first high-resolution X-ray spectrum of a putatively single 
Wolf-Rayet star. 400 ks observations of WR6 by the XMM-iVeutfon-telescope 
resulted in a superb quality high-resolution X-ray spectrum. Spectral analysis 
reveals that the X-rays originate far out in the stellar wind, more than 30 stellar 
radii from the photosphere, and thus outside the wind acceleration zone where the 
line-driving instability could create shocks. The X-ray emitting plasma reaches 
temperatures up to 50 MK, and is embedded within the un-shocked, "cool" stellar 
wind as revealed by characteristic spectral signatures. We detect a fluorescent Fe 



-2 - 



line at ~ 6.4 keV. The presence of fluorescence is consistent with a two-component 
medium, where the cool wind is permeated with the hot X-ray emitting plasma. 
The wind must have a very porous structure to allow the observed amount of 
X-rays to escape. We find that neither the line-driving instability nor any alter- 
native binary scenario can explain the data. We suggest a scenario where X-rays 
are produced when the fast wind rams into slow "sticky clumps" that resist ac- 
celeration. Our new data show that the X-rays in single WR-star are generated 
by some special mechanism different from the one operating in the O-star winds. 

Subject headings: Stars: winds, outflows — Stars: Wolf-Rayet — Stars: individ- 
ual: WR 6 — X-rays: stars 



1. Introduction 



Massive stars reach the Wolf-Rayet (WR) evolutionary phase when their hydrogen fuel 
has been consumed or lost and the products of nuclear fusion appear in their atmosphe res, be- 



fore e nding their short lives among the progenitors of core-collapse supernovae (SNe) (ISmartt 



20091 ). These stars are millions of times more luminous than the sun and drive strong stel- 
l ar winds by radiation pressure exerted through a bsorption and scattering in spectral lines 



(INugis Lamers 



2002 



Grafener fe Hamannll2005f) . This "line-driving" mechanism is known 



to be unstable (ILucy &: Whitdll980l ) and is thought to create hydrodynamic shocks that heat 
the plasma to a few million degrees with the emission of thermal X-rays (IFeldmeier et al. 



i_pi 

3). 



19971 ). For O-type stars, whic h have weaker winds than WR stars, this m odel is largely con 



sistent with the observations (jOskinova et al.l 120061 ; iGiidel fc Nazel 120091 ). The winds of WR 
stars are also expected to suffer from the line-driving instability (LDI) (IGayley fc Owocki 
19951 ) and thus may emit X-rays similar to the O-type stars. This conjecture was not verified 
by observations so far. 

In this Letter we report the first high-resolution X-ray spectrum of a putatively single 
WN-type star and its analysis. 



The WN star WR 6 



The object of our study, WR6 (EZ CMa, HD 50896), has been successfully modeled 
as a hydrogen-free WN s tar of subtype WN5 by fitting its UV, visible, and infrared spec- 
tra (jHamann et al.ll2006l ). Its stellar temperature is 90 kK, and its luminosity is 10 5 ' 6 L Q 
when adopting a distance of 1.82 kpc as implied from its membership to the association 
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Collinder 121. The mass-loss rate is M ~ 2 x lO _5 M yr _1 , and the wind reaches a ter- 
minal velocity of ~ 1,700 km s _1 . Hence the wind carries a mechanical power Mv 2 OQ j2 
equivalent to 1O 4 L , which is about 3% of the stellar luminosity. 



The WR winds are strongly inhomogeneous (ILepine fc Moffatl 119991 ) . Narrow spectral 
features drifting over the b road emission- line profiles are interpreted as showing radially 
accelerated blobs of matter (jMoffat et al.lll988l ). Quasi-periodic variability could be related 
to rotation ( iDessart fc Chesneaul 120021 ) and caused by spiral-like density patterns in the 
wind roote d in the photosphere, similar to the corotating interaction regions observed in the 
solar wind (lMullanlll984[ iMorel et al.lll997l ) . WR 6 shows considerable photometric, spectral, 
and p olarimetric variability on the time scale of 3.766 d (IDuijsens et al.lll996t ISt-Louis et al. 
20091 ). 



I ts thermal radio spec trum (IDougherty &: Williams! 120001 ) and complex variability pat- 



tern ( IDuijsens et al.l Il996h support the conclusion that WR6 is not a binary but a single 



star. Previous X-ray observations of WR6 at lower spectral re solution confirmed that its 



2003; 


Oskinova 


2005; 


Skinner et al. 


2010) 



3. Observations and Spectral Modeling 

The X-ray data on WR6 were taken with the X-Ray Multi-Mirror Satellite XMM- 
Newton. Its telescopes illuminate three different instruments which always operate simulta- 
neously: RGS is a Reflection Grating Spectrometer, achieving a spectral resolution of 0.07 A; 
RGS is not sensitive for wavelengths shorter than 5 A. The other focal instruments MOS and 
PN cover the shorter wavelengths; their spectral resolution is modest (E/AE «20 - 50). 

The data were obtained at four epochs in 2010 (Oct. 11 and 13, Nov. 4 and 6). The total 
exposure time of 400 kilosecond was split into four individual parts. Our data reduction 
involved standard procedures of the XMM-Newton Science Analysis System v. 10.0. 

The RGS spectrum is shown in Fig.[TJ The X-ray luminosity of WR6 in the 0.3 - 
12.0 keV band is L x ~ 8 x 10 32 ergs _1 , or O.1L . This constitutes 10~ 5 of the wind's 
mechanical power. The RGS spectrum is dominated by strong and broad emission lines 
of metals in accordance with the WN-wind abundances characteristic of CNO processed 
material. 



To model the observed X-ray spectrum we generated emissiy ities for a hydrogen- free 
plasma with the APEC thermal plasma model ( jSmith et al.l 120011 ). For most elements the 
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abundances were fixed to the values derived from the UV and optical spectra. 

In this stage we approximate the line profiles by Gaussians, with a common Doppler shift 
against the laboratory wavelength of — 650kms _1 and a width (FWHM) of 3000 km s -1 as 
the best values found by the automatic procedure from fitting the RGS data. The optim um 



fits to the observed spectra were obtained with the ISIS software (INoble fe Nowakl 120081 ) . 



The observed spectra, the best-fitting model, and the residuals are shown in Fig.|5J In 
order to achieve an acceptable fit, it was necessary to compose the emission from plasmas of 
three different temperatures (1.6, 7.0 and 45 MK) with respective emission measures of 253, 
86, and 22 x 10 54 cm -3 . The hottest component is required to reproduce the 2-5 A continuum 
and the emission lines from Fe xxv, Ca xix, Ar xvil, and S XV. Nevertheless, there are still 
residuals at some lines, such as Mg XII and Ne x, indicating that a three-temperature model 
is not fully sufficient, or that some abundances are not well determined. O and C lines are 
absent or very weak, as expected for a plasma with the chemical composition of a WN-type 
wind. 



To account for the absorption within the wind we applied the vphabs model ( lArnaud 



19961 ). modified for the absence of hydrogen. Strong absorption, significantly exceeding that 
of the interstellar medium towards the star, is evident. The ratio of fluxes in the Nvi and 
Nvn Lya lines is reproduced when absorption in the cool wind is included, especially the 
K-shell ionization edge of Niv that is located between the two aforementioned lines (see 
FigE}. 

The low-resolution EPIC spectra provide data up to 12keV. Significant residuals are 
encountered in the region around the Fe xxv line at about 1.9 A (see Fig. [2]). The observation 
shows a definitely broader feature than reproduced by the plasma model. Obviously, this 
complex includes a further blend. The fit of an additional gaussian profile yields a flux 
of 4.5 x 10~ 7 photons cm -2 s _1 and a line center at 1.92 A. The latter value agrees within 
the uncertainties with the wavelength of the Fe Ka line (1.94 A). Hence we conclude that 
cool-wind material must be irradiated by X-rays which are hard enough (> 7.1 keV) to 
excite this line by fluorescence. Our PoWR wind models show that in the outer wind the 
leading ionization stage of iron is FeV. The presence of fluorescence is consistent with a 
two-component medium, where the cool wind is permeated with the hot X-ray emitting 
plasma. 
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The cool wind model of WR 6 



To determine the physical conditions in the "cool" component of t he stellar wind, we 
employ the PoWR model atmosphere code (IHamann fc Grafenerl 12004 ) . The PoWR code 
solves the non-LTE radiative transfer simultaneously with the equations of statistical and 
radiative equilibrium. Complex model atoms with thousands of transitions are taken into 
account. The extensive inclusion of the iron gr oup elements is impo rtant because of their 
blanketing effect on the atmospheric structure ( jGrafener et al.ll2002l ). A particular stellar 
atmosphere model is defined by the effective temperature, surface gravity, luminosity, mass- 
loss rate, wind terminal velocity, and chemical composition. 

For the supersonic part of the wind, we adopt the radial dependence of velocity as 
v(r) = Voo(l — 1/r), where r is in units of the stellar radius i?*, and the terminal velocity 
Voo is a free parameter. In the subsonic region, the velocity field is defined such that the 
hydrostatic density stratification is approached. The PoWR models account for stellar- wind 



clumping in the standard yo lume-filling factor 'microclumping' approximation (lHilrierlll991 
Hamann fc Koesterkdll998l ). Our best model of WR6 has a clumping factor D = 20, where 
D = (p 2 )/(p) 2 . The synthetic spectra are calculated over the whole spectral range from UV 
to IR, and then compared to the observed spectra and photometric fluxes. 

With the final wind model being adopted, the stratifications of the density, opacity, 
ionization, and radiative flux are defined. The mass absorption increases with distance from 
the star, partly because the fully ionized helium recombines to He II with its strong bound- 
free opacity. This increase of the mass absorption coefficient partly compensates for the 
decrease of matter density, keeping the wind opaque to large radii. 

In Fig.[3] (upper panel) we show the radius in the wind where the continuum optical 
depth on the radial ray becomes unity, as predicted by PoWR model. Roughly speaking, 
only X-rays emitted from outside that radius can escape from the wind and be seen by a 
distant observer, unless inhomogeneities create some porosity ("macro clumping"). 



5. The analysis of line ratios in He-like ions 



Helium-like ions show a group of three lines, consisting of a forbidden (/), an intercom- 
bination (i), and a resonance (r) transition - the so-called fir triplet. The fir triplets of 
Si xill, Mg XI, Ne IX, and N vi are present in the RGS spectrum of WR6. As illustrated 
in Fig.Hl the line ratios observed in WR6 are very different from those observed in typical 
single O-type stars. For N vi the ratio fji exceeds un ity in WR6, while O stars show f ji < 1 
(ILeutenegger et al.ll2006l ; IWaldron fc Cassinellil 120071 ). 
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For each He-like ion, the ratio of fluxes between the forbidden and the intercombination 



comp onent, R, is sensitive to the electron density and the ultraviolet flux ( IBlumenthal et al. 
1972h : 

R(r) = — (1) 

[) l + 0(r)/0 c + iV c (r)/iV c ' {) 

is the photoexcitation rate from the term 2s 3 S to 2p 3 P, and N e is the electron 



where 

density. The quantities IZn, 



and N p depend only on at omic parameters and the electron 



temperature (IBlumenthal et al.lll972t iPorquet et al.ll200ll ). 

We use the PoWR model to calculate the values of R(r) for Mg XI, Ne IX, and N VI 
as function of the radial location in the wind of WR6. The density in the X-ray emitting 
shocks is not known, therefore we computed two sets of R(r), one neglecting the density 
term, and another one assuming that the density in the shock is the same as in the ambient 
"cool" wind but the plasma is fully ionized. According to the s trong-shock condition, th e 
hot gas density would be even a factor of four higher than that (jZeldovich &: Raizerlll966l ). 
The photo-excitation rates 4>(r) are computed at each radius from the radiation intensity as 
provided by the PoWR model, which accounts not only simply for geometric dilution, but 
also for the diffuse radiative field. Fig. [3] (low panel) shows for the N VI triplet the predicted 
R ratio as function of the radial location of the emitting plasma with the best-fit value of 
~ 200 R* when the collisions are accounted for. The N vi emission line indicates a plasma 
temperature of ~ 1.58 MK, which could be reached in a strong shock with a velocity jump 
by ps 400 km s -1 . 



6. Modeling of X-ray emission line profiles 



The broadening of the X-ray lines from WR6 is consistent with the known terminal 
velocity of the wind of ~ 1700kms _1 . The lines centers are displaced from their laboratory 
wavelengths by ~ —0.06 A (see Fig. [5]). Such blue-s hifts are predicted if th e line emissio n 
is distributed within a partly absorbing wind (e.g. iMacfarlane et al.l Il99ll ; llgnacd l200fl ). 
However, according to the standard model, the wind of WR6 remains optically thick till 
large distances from the star (550 - 900 i?* for A > 20 A). X-rays can hardly emerge from 
interior to this location unless the wind is very porous. 



We employed our 2-D stochastic wind code (jOskinova et al.ll2004l ) to model the observed 
X-ray lines. Relaxing the microclumping approximation, the clumps may have arbitrary 
optical depth ("macro clumping"). This wind fragmentation alters the r adiative transfer 



drast ically, compared to a homogeneous wind with the same mass-loss rate (IFeldmeier et al. 



20031 ). The X-rays can escape from deeper in the wind than indicated in the upper panel in 
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Fig.[3]for the smooth- wind model. 

We assume that parcels of X-ray emitting gas are distributed between clumps of continuum- 
absorbing, cool wind. The wind is fully fragmented throughout the whole relevant radial 
range. The X-ray emission is placed only at radii larger than r em . Since thermal X-ray 
emission arises from the decay of collisionally excited levels and from free-free emission, is 
plausible to assume that the emissivity scales as the square of the density, r]\ oc p 2 . With this 
scaling, most X-rays are released close to r em . The intensity along any ray is reduced by the 
factor exp(— r£), where r£ is the summation over the optical depths through all individual 
clumps which are encountered. 

A grid of line models was calculated for different values of the radius r em and (N c ) - 
the time-averaged number of fragments in a radial direction. The radial dependence of the 
mass-absorption coefficient and the wind opacity was accounted for in the models. Even with 
porosity, emission arising from too close to the star is totally trapped, while the adopted 
scaling of the emissivity with p 2 requires an unrealistic amount of energy being converted 
into X-rays. This excludes small values of r em . In the other extreme, when X-ray emission 
starts only at radii where the wind is already nearly transparent, the observed blue-shift and 
asymmetry of the line profiles cannot be reproduced. The model line profiles shown in Fig. [5] 
were calculated with r em = 30/2* and (N c ) = 50. 



n The Origin of X-ray Emission in WR stars 



7. O 



Observed temperatures up to 50 MK require shock velocity jumps comparable to the 
wind speed. From UV spectra of WR6 there is direct evidence that a fraction of the 
wind flows wit h much higher ve locity (up to 3,100 km s _1 ) than the general wind speed 
(1,700 km s _1 ) (IPrinja et al.lll990l ). We suggest that the hottest gas is produced when very 
fast flows impacts on slowly moving clumps within the wind acceleration zone. Such extreme 
shocks cannot be maintained to large radii, but the hardest X-rays can escape from the inner 
wind owing to smaller absorption cross-sections at higher energies. 

Exceptionally large, pancake-shaped clumps which form close to the photosphere at 
low velocities would resist the radiative acceleration. Being an obstacle to the generally 
much faster wind, they create strong reverse shocks that may add more mass to the clump 
if the gas cools quickly, or generate bow shocks if it does not (ICassinelli et al.l 120081 ) . Self- 
consistent kinematic models of Gayley (in prep.) show that 500kms _1 shocks may be present 
at ~ 50 i£* in the wind, suitable for X-ray escape if the wind is sufficiently porous. 
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8. On the binary hypothesis for WR6 



The X-ray emission of WR6 varies considerably (lPollock!ll987l ; IWillis fe Stevens! Il996l ). 
Our new observations confirm and extend these results. The XMM-Newton light curve shows 
different variability time scales, including a variation of 20% over a month's span. We also 
observe some spectral changes. The X-ray variab i lity d oes not follow the 3.766 d period 
which is seen in the optical and UV ( 1 Willis et al.l Il989l ) . The X-ray variability of WR 6 
further highlights the differences between single WR and O stars. The latter are remarkably 
constant X-ray sources (Naze et al. 2012, submitted). 



The optical variability of WR6 has long been known ( )Wilsonlll948D, and it appears t o 



be incompatible with stellar pulsations or the Be star phenomenology flRobert et al 



Most plausibly, th e variability is due to changes of wind structures (jSt .-Louis et al. 



1992 . 



1993 



Flores et al.l 1201 if ). This explanation is in line with our new X-ray observations. 



Firmani et al.l ( Il980l ) established that if the 3.766 day period were associated with bi- 
narity, the mass of the secondary star could not be larger than 1.5 M & . A star of such 
mass and age (assuming coeval formation of both binary components) would be either a 
degenerate neutron star (NS) or a pre main-sequence young stellar object (TTauri type 
star). 

The NS scenario was criticized because of t he relatively low X-ray luminosity of WR 6 
and the lack of evidence in its UV spectrum ( IStevens fc Willisl Il988[ ). The WR+TTau 
scenario is also extremely unlikely, as the high- resolution X-ray spectrum of W R6 is incom- 
patible with the X-ray properties of young stellar objects (IKastner et al.ll2002l ). 



Variability is a common property of WR stars. Among rigorousl y monitored WN stars, 
40% show optical variability similar to WR6 ( jChene fc St-Louisll201ll ). The assumption that 
40% of all WN stars have a NS or a T Tauri type companion is not realistic. 

Our high-resolution X-ray spectra show that the X-ray emitting plasma moves at about 
the same velocity as the cool wind. The X-ray line blue-shifts do apparently not change with 
time. The X-ray emission line spectrum is compatible with the WN star abundances. All 
these facts add further evidence against WR6 being a binary system. 



9. Conclusions 

Assuming that WR 6 is typical, our results show that in the winds of WN stars strong 
shocks are active at large radii from stellar photosphere. A possible explanation involves slow, 
dense clumps that accumulate from the wind close above the photosphere without being ef- 
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fectively accelerat ed. The seeding of such slow cl umps might be associated with (sub) surface 



convective zones (jCantiello fe Braithwaitd l201ll ). The effects of such wind anisotropics on 



the distribution of circumstellar matter could influence the behavior of SN Ibc. 

In any case, from high-resolution X-ray spectral data of WR 6 we conclude that unknown 
mechanisms must operate in WR winds. Identifying these mechanisms poses a challenging 
problem for the theory of stellar winds. 



Based on observations made with XMM-Newton, an ESA science mission with instru- 
ments and contributions directly funded by ESA member states and the USA (NASA). We 
thank R.K. Smith and A. Foster for help in computing emissivities for a hydrogen-depleted 
plasma. DPH was supported by NASA through the Smithsonian Astrophysical Observatory 
contract SV3-73016 for the Chandra X-Ray Center and Science Instruments. Funding for 
this research has been provided by DLR grant 50 OR 1101 (LMO). 
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Fig. 1. — The RGS spectrum of WR6 with the strongest emission lines identified. The error 
bars correspond 3cr. The RGS 1 spectrum is shown from 18.2 A to 26.7 A, otherwise the 
RGS 2 spectrum is plotted 
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Fig. 2. — The X-ray spectrum of WR6, obtained with different detectors: RGS (top panel), 
MOS (middle), and PN (bottom). Black curves display the observations; the model is shown 
in red. The lower part of each panel displays the residuals. The PN panel includes an inset to 
the right showing the transmission of the cool wind plus interstellar absorption for A > 20 A. 
The N IV K-shell absorption edge is located between the N vi and N VII lines (red vertical 
bars). An insert to the left in the PN panel shows a broad spectral feature that is not 
reproduced by the thermal plasma model (blue dotted line). An additional line at 1.92 A, 
which we identity as fluorescent FeKa, improves the fit (red histogram). 
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Fig. 3. — Upper panel: radius where the continuum optical depth reaches unity, as predicted 
by the PoWR model for the "cool" wind component of WR6, when inhomogeneities are only 
accounted for in the "microclumping" approximation. Lower panel: Dependence of the line 
ratio R — f/i for the N vi lines as a function of the radial location of the emitting plasma. 
Based on the PoWR model, the red curve accounts only for radiative de-population, while 
the blue curve includes also collisional process under the assumption that the shocked plasma 
has the same density as the cool wind. The measured value is shown as horizontal green 
line, with the green shaded area representing the 3a confidence band of the measurement. 
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Fig. 4. — Upper panel: RGS spectrum of WR6 in the range of the Nvi triplet of lines 
(black). The vertical error bars correspond to 3cr. The vertical lines indicate the restframe 
wavelength of the resonance (A r ), intercombination (A;), and forbidden (Af) line, respectively. 
The solid red line shows a fit to these lines by a formal model with three Gaussian profiles. 
Lower panel: Same as upper, but for the 041 star ( Puppis. The red line shows the scaled 
Gaussian fit from the upper panel. 



- 16 - 



3200 h-NVII 



2400 - 



§ 1600 

o 
o 



800 - 



l — i — i — r 



j L 




24.5 



24.6 



24.7 



24.8 



24.9 




14.8 



14.9 15.0 

Wavelength (A) 



15.1 



Fig. 5. — Upper panel: Part of the spectrum of WR6 centered on the N VII Lya line. The 
observed spectrum is shown as a histogram (black). The laboratory wavelength (Ao) of 
N vilLya A 24. 779 A is indicated by a vertical dashed line. The model line (smooth red) is 
a blend of N vilLya A 24. 779 A and A 24. 785 A with a ratio of line emissivities 1 : 0.5. Lower 
panel: The same as upper panel, but for Fexvil A15.014A. The laboratory wavelength 
(Ao = 15.014 A) is indicated by a dashed line. The model line (smooth red) is a blend of 
Fexvil A 15.014 A and A 14.892 A with a ratio of line emissivities 1 : 0.1. 



